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Abstract 

In this paper, we generalize our results in [5] to triangulated surfaces in hyperbolic 
background geometry, which means that all triangles can be embedded in the stan¬ 
dard hyperbolic space. We introduce a new discrete Gaussian curvature by dividing 
the classical discrete Gauss curvature by an area element, which could be taken as the 
area of the hyperbolic disk packed at each vertex. We prove that the corresponding 
discrete Ricci flow converges if and only if there exists a circle packing metric with zero 
curvature. We also prove that the flow converges if the initial curvatures are all nega¬ 
tive. Note that, this result does not require the existence of zero curvature metric or 
Thurston’s combinatorial-topological condition. We further generalize the definition 
of combinatorial curvature to any given area element and prove the equivalence be¬ 
tween the existence of zero curvature metric and the convergence of the corresponding 
flow. 


1 Introduction 

Consider a compact surface X with a triangulation T = {G, i?, F}, where the symbols 
V, E, F represent the set of vertices, edges and faces respectively. A positive function 
r : y —>■ (0, -boo) defined on the vertices is called a circle packing metric and a function 
^ : E ^ [0,7r/2] is called a weight on the triangulation. Throughout this paper, a 
function defined on vertices is regarded as a column vector and N = is used to denote 
the number of vertices. Moreover, all vertices, marked by vi,...,V]\f, are supposed to be 
ordered one by one and we often write i instead of u*. Thus we may think of circle packing 
metrics as points in ^ times of Cartesian product of (0, oo). 

Given (A, T, $), every circle packing metric r determines a piecewise linear metric on 
X by attaching edge etj a length bj = -b -b 2rirjCos{^ij). This length structure 
makes each triangle in T isometric to a Euclidean triangle. Furthermore, the triangulated 
surface (A, T) is composed of many Euclidean triangles glued coherently. This case is 
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called Euclidean background geometry in [T]. If the length lij of the edge is determined 
by the hyperbolic cosine law cosh lij = cosh r* cosh rj + sinh sinh rj cos ^ij , then each 
triangle in T is isometric to a hyperbolic triangle embedded in the standard hyperbolic 
space In this case the triangulated surface {X,T) is composed of many hyperbolic 
triangles glued coherently, which is called hyperbolic background geometry. 

We use {X,T, <I>,1H^) to denote the space we want to study in the following, where X 
is a closed surface, T is a fixed triangulation of X, is a fixed weight function defined on 
edges, and represents the hyperbolic background geometry. It was proved by Thurston 
[7] that, whenever {i,j,k} G F, these three positive numbers lij,lik:ljk satisfy the trian¬ 
gle inequalities. Thus the combinatorial triangle {i,j,k} with lengths lij,lik,ljk forms a 
hyperbolic triangle in H^. 

For the hyperbolic triangle AviVjVk, the inner angle of this triangle at Vi is denoted 
by and the classical combinatorial Gauss curvature Ki at Vi is defined to be 

K, = 27r- ( 1 - 1 ) 

Notice that can be calculated by hyperbolic cosine law, thus and Ki are elementary 
functions of the circle packing metric r. 

Obviously, discrete Gaussian curvatures K is determined by circle packing metric r. 
Gonversely, Thurston proved [7] that ii'-curvature map r ^ K \s injective, which implies 
that the metric r is determined by its X-curvature. Among all circle packing metrics, a 
metric with all curvatures zero, which is unique if it exists, is of special interest. Zero 
curvature metric does not always exist. At first glance, the existence of zero curvature 
metric depends on the topological information of the surface. In fact, for the hyperbolic 
triangle AviVjVk, we have -|- 0*^ -|- 0]^ = vr — Area{AviVjVk)- Using this formula, one 
get the combinatorial Gauss-Bonnet formula [T] 

N 

Y,Ki = 27rxiX)+Area{X), ( 1 . 2 ) 

i=l 

which implies, if there exists a circle packing metric with zero curvature, then x{Ai) must 
be negative. However, only topological constraint is not enough. In his work on construct¬ 
ing hyperbolic metrics on 3-manifolds, Thurston [7j discovered additional combinatorial- 
topological conditions which is equivalent to the existence of zero curvature metric. 

Proposition 1.1. ([7]) Given (A, T, <h,H^), there exists a zero curvature metric if and 
only if the following two combinatorial-topological conditions are satisfied simultaneously: 

1. For any three edges ei, 62 , 63 forming a null homotopic loop in A, if ^ '^1 

then ei, 62,63 form the boundary of a triangle of T; 
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2. For any four edges ei, 62 , 63 , 64 forming a null homotopic loop in X, if ‘h(ej) > 
27r, then 61 , 62 , 63,64 form the boundary of the union of two adjacent triangles. 


Bennett Chow and Feng Luo introduced the combinatorial curvature flow to study the 
problem. They defined a combinatorial Ricci flow 


-^ =-RTj sinh rj (1.3) 

in [T] and proved that the flow (II.3p converges if and only if there exists a metric with 
all curvatures are zero, if and only if Thurston’s combinatorial-topological conditions are 
satisfied. 

Inspired by [T] and [ 2 ], we studied the negative gradient flow of discrete Calabi energy. 


—d- = — ARTjsinhrj, 
dt 


(1.4) 


which is called combinatorial Calabi flow in [ 2 ]. We proved that combinatorial Calabi flow 
converges if and only if zero curvature metric exists, assuming a uniform upper bound of 
metric r. 

The paper is organized as follows. In Section O we introduce a new definition of 
discrete Gauss curvature for triangulated surfaces in hyperbolic background geometry and 
the corresponding combinatorial Ricci flow. In Section [31 we prove that the existence 
of zero curvature metric is equivalent to the convergence of combinatorial Ricci flow we 
introduce in the paper, which provides another description of Thurston’s existence of circle 
packing. In Section we prove that the existence of metrics with nonpositive curvatures 
ensures the existence of zero curvature metric. In Section [5l we generalized the definition of 
combinatorial curvature to general area element and then prove that the existence of zero 
curvature metric in this sense is also equivalent to the convergence of the corresponding 
combinatorial curvature flow. 


2 Definitions of discrete Gaussian curvature and Ricci flow 


By the analysis in [5], the classical combinatorial Gaussian curvature is not a good can¬ 
didate for approximating the classical Gaussian curvature. Following the idea and the 
approach in [5], we have the following definition of combinatorial Gaussian curvature. 


Definition 2.1. Givin (A, T, <I>,H^) with circle packing metric r, the modified combina¬ 
torial Gaussian curvature at the vertex i is defined as 


Ri = 


Kj 

dvr sinh^ ^ ’ 


( 2 . 1 ) 
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where Ki is the classical combinatorial Gaussian curvature given by dLH). 

In fact, the denominator dvrsinh^ ^ is just the area of the hyperbolic disk with radius r 
packed at vertex i. To study the constant combinatorial curvature problem, we introduce 
a combinatorial Ricci flow which is an analogue of the smooth Ricci flow. 


Definition 2.2. Given (X, T, d>,lH^) with circle packing metric r, the combinatorial Ricci 
flow is defined as 

= -R^9i, ( 2 . 2 ) 

where gi = sinh^ 


For the following application, we usually write the combinatorial Ricci flow (I2.2p in 
the following form 


dui 

dt 



(2.3) 


where Ri = -—— is a modification of the curvature Ri and m = In tanh ^ is a coordi- 
nate transformation, which maps R^q homeomorphically onto R^q- 

The flow (j2.2p is in fact an ODE system. Therefore, the solution always exists locally 
around the initial time t = 0 since all coefficients are smooth and then locally Lipschitz 
continuous. It’s easy to get 


Proposition 2.3. If the hyperbolic combinatorial Ricci flow (j2.2p converges, then there 
exists a circle packing metric with 0 R-curvature. 


Proof. Suppose the solution r{t) —)■ r* as t —>■ +oo. Then we have g{t) —>■ g*, which 
implies that there exists a sequence —)■ +oo such that 


gi{n + 1 ) - gi{n) = g[{in) = -RmiCn) 0 


as n —>■ + 00 . As r{t) r*, we have R{t) —>• R*. Thus R*g* = 0, which implies R* = 0 
and r* is a metric with zero curvature. □ 


Obviously, the X-curvature, the R-curvature and the R-curvature are different. How¬ 
ever, a metric has zero X-curvature if and only if it has zero R-curvature, if and only if 
it has zero R-curvature. Thus in the following we needn’t distinguish them and just call 
them “zero curvature metric” for short. 

To study the long time behaviors of combinatorial Ricci flow (12.211 . we shall use two 
different ways in the following subsections. 
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3 Convergence of Ricci flow when zero curvature metric 
exists 


In this subsection, we shall study carefully the relationship between the existence of zero 
curvature metric and the long time behavior of the combinatorial Ricci flow. We also 
derive some convergence results under the assumption of existence of a zero curvature 
metric, which is necessary for the convergence of flow (12.2p by Proposition 12.81 The main 
tools used here are standard in ODE theory. 

The following two Lemmas are useful. 


Lemma 3.1. m [ 3 ]) Given {X, T, $, H^) with circle packing metric r. u is the coordinate 
transformation of r with m = In tanh ^ ■ Then 


,Kn) 

d{ui, ■ ■ ■ ,un) 


A + Lb 


is positive definite, where A is a positive diagonal matrix and Lb is positive semi-definite 
with rank N — 1 and kernel t(l, • • • , 1)^ and could be written as 




k'^i 

-Bij , j ~ f, 
0 , 


with Bij > 0. 


□ 


Lemma 3.2. (Lemma 3.5, [T]) For any e > 0, there exists a number I so that when > I, 
the inner angle 9i in the hyperbolic triangle AviVjVk is smaller than e. 

Proof. This result has been stated in [T], here we give a proof just for completeness. 
We will prove that 0 uniformly as r* —>■ -|-oo. By the hyperbolic cosine law, we have 

jh cosh lii cosh lii^ — cosh 

- - 

Siriil l-ij SlRll (jik 

_ cosh{lij + Ijk) + cosh{lij - Ijk) - 2cosh/jfc 

cosh(/jj + lik) - cosh(/jj - kk) 

1-A ’ 


where A = kA ^ coshijk ^ ^ ^ 

cosh{lij+lik) ^ 


0 < A < 


cosh{lij+lik) ■ 
cosh lij 


< 


1 


< 


1 


cosh(^jj -|-cosh lik coshr. 
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we have A —>■ 0 uniformly as r* — >■ +cxd. So we just need to prove that ^ ^ 0 uniformly as 
Vi —>■ + 00 . Note that 

coshljk < coshrj coshr^ + sinhrj sinbr^ = cosh(rj + r^), 
we just need to prove that lij + Uk — {vj + r^) ^ +oo uniformly as ^ +oo. Note that 

cosh lij = cosh Tj cosh Vj + sinh rj sinh rj cos > cosh r* cosh rj > — , 


we have 


hj > In cosh lij > Tj + rj — In 4. 


Similarly, we have 


kk > In cosh lik >ri + rk-In 4. 


Then we have 

kj + kk - (rj + rfc) > 2rj - 2 In4 +oo 


uniformly as rj —?> +oo. 


□ 


Theorem 3.3. Assuming there exists a zero curvature metric r* on (X, T, Then 

the solution r{t) to the combinatorial Ricci flow (I2.2h exists for all time t € [0, +oo) and 
converges exponentially fast to r*. 


Proof. Denote u* as the corresponding u-coordinate of r*. Consider the combinatorial 
Ricci potential (first introduced by Chow and Luo in [1]), 


F{u) 



(3.2) 


Then Hess{F) = L is positive definite. We had proved that F{u) > F{u*) = 0 and u* is 
the unique minimum point and zero point of F. Moreover, lim F{u) = Too (Lemma 

|| u ||—>-+00 


B.l, [3]). Then 


dF{u{t)) 

dt dt 




E 


Kf 


27r sinh^ n 


< 0 . 


(3.3) 


So F is decreasing along the flow. Using lim F{u) = +oo, we know {u(t)} is bounded 

|| tt ||—>-+00 

in R< 0 ) i-C- there exists c > 0, such that Ui{t) > —c for any i £ V and t G [0,+oo). 
Therefore rj(f) > > 0 for any i £ V and t £ [0, +oo), where 5 = In ■ 
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The combinatorial Ricci flow (12.21) can be written as (coshrj)'(t) = —Note that Ki 
is uniformly bounded by (2 — d)Tr < Ki < 2 tt, where d is the maximal degree at vertices. 
Hence coshrj(t) < coshrj(O) + and then rj(t) < ln2 + In (coshri(O) + for each 

i €V, which implies that the solution r{t) exists for all time t G [0, +oo). 

We claim that r(t) is uniformly bounded above. If not, then there exists at least one 
vertex i £V, such that lim rRt) = +oo. For vertex i, using Lemma 13.21 we can choose 

t—>-+oo 

I > 0 large enough so that, whenever r* > I, the inner angle 0^^ is smaller than where 
di is the degree at vertex i. Thus Ki = 2^ — > vr. Choose a time to such that 

f'ii'to) > ^ this can be done since lim rj(t) = +oo. Denote a = inf{ t < to I ri{t) > I }, 

t^+oo 

then rj(o) = 1. Let’s look at what happens to flow ()2.2p in the interval [a, to]. Note that 
flow (|2.2p can be written as = — 2 -K^nhr- ’ ~Ki < —vr when a < t < to, then 
r'(t) < 0 and hence rj(t) < ri{a) = I, which contradicts rj(to) > L 

The arguments above show that {r(t)|t G [0, +C)o)} is compactly contained in R^o- 
show there exists a sequence t Too, such that r(^„) —> r*. In fact, F(u(+oo)) exists. 
Furthermore, there exists a sequence t Too such that 

F(n(n + l))-FKn)) = (T(u(t)))'|^„ = = -J] 27 rsinh 2 r,.R 2 |^^ ^0^ 


Since {r(t)} CC R> 0 ) ^(^n) have a convergent subsequence which is still denoted as r(^„) 
with r{^n) —t T^*- Then we get R{r*) = 0, and hence K{r*) = 0 and R{r*) = 0. By 
Thurston’s injective property of iL-curvature map, r* is the unique zero curvature metric. 
By calculating the differential of the right hand of (12.31) at u*, we get 


DI 




where S = dm^ilsinh^ ri, • • • ,sinh^rAr}. As —= —S 2 S“ 2 LS “2 ~ S” 

D\u*{—R) has N negative eigenvalues, which implies that u* is a local attractor of (12.31) . 
Combining this with the fact that u{^n) converges to u*, we get u{t) —)• u* for any initial 
metric u(0) with exponential convergence rate by ODE theory. □ 


1 

2 ~ 


■5LS- 


4 Convergence with nonpositive initial curvature metrics 

In the last subsection we proved combinatorial Ricci flow converges exponentially fast to 
zero curvature metric if there exists a zero curvature metric. It’s interesting that we can 
get a convergence result without the assumption of the existence of zero curvature metric. 
The main tool used in this subsection is a discrete version of maximum principle first 
developed by the authors in [5]. 
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We first derive the evolution of modified curvature Ri along the combinatorial Ricci 
flow. 

Lemma 4.1. Along the combinatorial Ricci flow (|2.2I) . the modihed combinatorial curva¬ 
ture Ri evolves according to 


dRi 

dt 


1 


{LR)i + 2 cosh TiR^. 


27r sinh^ r 

Proof. First note that Ui = Intanh we have 

d . , d 

— = smhri —. 

OUi OTi 

Then we have 

dRi I dKi K, . ^ d . . 1 dKi 

duj 2 tt sinh ri duj 2 tt sinh ri drj 2 tt sinh r* duj 

which implies that 
dRi 

dt 


(4.1) 


— 2Ri sinhrj coth ridij, 


^ > dRi duj 


duj dt 27r sinh^ rj duj 


^ IZ|l“^i + 2coshrjRf. 


□ 

To get the estimation of curvature Rj, we can use the following discrete version of 
maximal principle introduced in [5]. 

Proposition 4.2. (Maximum Principle) Let / : R x [0, T) ^ R be a function such 
that 

dfi 


dt 


> A/j + chj(/j), V(z,t) gRx [0,r) 


where A/j = ~ fi) with aij > 0 and 4>j : R —R is a local Lipschitz function. 

Suppose there exists Ci G R such that /j(0) > Ci for all i G R. Let if be the solution to 
the associated ODE 

' = ^,{ip) 


then 


dt 

¥^(0) = Cl, 


fi{t) > if{t) 


for all (z, t) £ V X [0, T) such that ip{t) exists. 

Similarly, suppose f : V x [0, T) —)• R be a C^ function such that 


dfi 

dt 


< A/, + chj(/j), V(z,t) gRx [0,r). 
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Suppose there exists (^2 G K, such that /j(0) < C 2 for all i G V. Let ip be the solution to 
the associated ODE 

I S = -!-<») 

\ m = C2, 

then 

fi(t) < p>(t) 

for all (i, t) € V X [0, T) such that V’(i) exists. 

Using the Maximum principle, we have the following corollary. 

Corollary 4.3. If Ri(0) < 0 for all i G V, then Ri{t) < 0 for all i € U along the 
combinatorial Ricci flow (12.21) . If Ri{0) > 0 for all i G V, then Ri{t) > 0 for all z G U 
along the combinatorial Ricci flow (|2.2p . 

Proof. As Ri and Ri always have the same sign, we just need to prove the corollary 
for Ri. As L = A + Lb hy Lemma I3.I1 we can rewrite the evolution equation for R as 

= ARi - —Ri + 2cosh.riR^, 

dt 27rsinh2ri 

where Ai?j = Yljr..^i^ij{Rj~Ri) with Ojj = 27rsinh^ r- ^ ^Pplyfog ^^6 maximum principle 
gives the proof of the corollary. □ 

Theorem 4.4. Given (A, T, <h, H^). If the initial curvature are all nonpositive, i.e. 
Ri(0) < 0 for each i G V, then there exists a zero curvature metric. Furthermore, the so¬ 
lution to combinatorial Ricci flow (j2.2p exists for all time and converges to zero curvature 
metric exponentially fast. 

Proof. Let r{t) be the solution. The uniform upper bound of r(t) follows the proof 
of Theorem 1,3.81 For the lower bound, note that the negativity of the modified curvature 
Ri{t) is preserved along the flow. Thus we have ^ = —Ri siuhr* > 0, which implies that 
u(i) > u(0) > 0. Hence we get {r(t)} CC H>o- Using Lemma [4.51 below, we get the 
proof. □ 

Lemma 4.5. If the solution r{t) to the combinatorial Ricci flow (12.2p stays in a compact 
subset of 1R.>0! then there exists a unique zero curvature metric. Moreover, r{t) converges 
exponentially fast to zero curvature metric. 

Proof. The condition {r(t)} CC implies the long time existence of the combina¬ 
torial Ricci flow. Consider the modified potential 
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where uq is an arbitrary fixed point in Then along the Ricci flow 


we have 




dt 


dt 


E 


K} 


2tt sinh r; 


2 — 


< 0 . 


(4.3) 


So G is decreasing along the flow (j2.2p . Generally, G is not bounded from below. However, 
G{u{t)) is bounded since {u(t)\t € [0,+oo)} is compact. Hence G(u(+oo)) exists. Then 
there exists a sequence —>■ +oo such that 

G{u{n + 1)) - G{u{n)) = {G{u{t))y\^^ = VG- = - ^ 27r.R^ sinh^ ^0. 


As ri(t) is bounded from above and away from zero, we have Ri{^n) —^ 0. Up to a sub¬ 
sequence, we can suppose that r(^„) —?> r* which has zero curvature. We can finish the 
proof by showing that D\u*{—R) = —has N negative eigenvalues and hence u* is 
a local attractor of flow (j2.3p or by using the conclusion of Theorem 13.31 directly. □ 


5 Combinatorial Ricci flow for general “area element” 

For arbitrary N positive functions Ai : R^q —)• (0,+oo), r Ai{r), we may take Ai as 
the most general form of “area element” at vertex i, where 1 < i < A^. We shall study 
combinatorial Gaussian curvature and combinatorial Ricci flow for this generalized “area 
element” Ai in the following. 

Definition 5.1. (A-curvature) The modified combinatorial Gaussian curvature with re¬ 
spect to area element Ai is defined as 

= (5.1) 

In the following we call it A-curvature for short. 


Obviously, a metric with zero A-curvature is exactly a metric with zero iF-curvature. 
Thus by mm, the zero A-curvature metric is unique if it exists. Inspired by the form 
(12.3p of combinatorial Ricci flow, we introduce the following A-flow to study A-curvature. 


Definition 5.2. (A-flow) The hyperbolic combinatorial Ricci flow with respect to area 
element Aj is dehned as 


dui 

dt 


-Ri, 


(5.2) 


where Ui = Intanh In the following we call it A-curvature flow or A-flow for short. 
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Theorem 5.3. Given (X, T, H^), ^-flow (15.21) converges if and only if there exists a 
zero curvature metric. Furthermore, if the solution of (15.21) converges, then it converges 
exponentially fast to the zero curvature metric. 

Proof. For the “only if’ part, the proof is the same as that of Proposition 12.31 For 
the “if’ part, the methods is the same as in the proof of Theorem 13.31 We only list some 
differences here. In the proof of Theorem 13.31 we first proved that all ri{t) are uniformly 
bounded below by a positive constant. This procedure can be used here again to get a 
uniform lower bound of r(t) as Ai is always positive. 

Suppose T is the maximal existing time of r{t) with 0 < T < +oo. We claim that 
ri{t) is uniformly bounded above for alH E F and t E [0, T). If not, there exists i E F 
and tn t T, such that ri{tn) t +oo. Choose the same I as in the proof of Theorem 13.31 and 
a no such that ri{tno) > I- Still denote a = inf{t < > 0) then ri{a) = 1. Note 

that A-flow satisfies r'(t) = —if*sinhrj/^j < 0, then rj(t) < rj(o) = I, which contradicts 
f'ii'tno) > I- Thus ri{t) is uniformly bounded above for all i and t, which implies that 
T = +00 and {r{t)\t E [0, +oo)} CC lR.>o- finish the proof by using Lemma 

which is still valid when we substitute sinhr^ by A^. □ 

Acknowledgements 

The authors would like to thank Dr. Wenshuai Jiang for reading the paper carefully 
and giving some valuable suggestions to improve the writing of the paper. The first 
author would also like to give special thanks to Dr. Yurong Yu for her supports and 
encouragements during the work. The research of the second author is partially supported 
by National Natural Science Foundation of China under grant no. 11301402 and 11301399. 
He would also like to thank Professor Guofang Wang for invitation to the Institute of 
Mathematics of the University of Freiburg and for his encouragement and many useful 
conversations during the work. 

References 

[1] B. Chow, F. Luo, Combinatorial Ricci flows on surfaces, J. Differential Geom, Volume 63, 
No. 1 (2003), 97-129. 

[2] H. Ge, Combinatorial Calabi flows on surfaces, Preprint at arXiv:1204.2930 [math.DG]. 

[3] H. Ge, X. Xu, 2-Dimensional Combinatorial Calabi Flow in Hyperbolic Background Ceometry, 
Preprint at arXiv:1301.6505 [math.DG]. 

[4] H. Ge, X. Xu, Discrete Quasi-Einstein Metrics and Combinatorial Curvature Flows in 3- 
Dimension, Preprint at arXiv:1301.3398 [math.DG]. 

[5] H. Ge, X. Xu, A combinatorial Yamabe problem on two and three dimensional manifolds. 
Preprint at arXiv:1504.05814 [math.DG]. 


11 



[6] F. Luo, Combinatorial Yamabe flow on surfaces. Commun. Contemp. Math. 6 (2004), no. 5, 
765-780. 

[7] Thurston, William: Geometry and topology of 3-manifolds, Princeton lecture notes 1976, 
http: / / www.msri.org/publications/books/gtSm. 

(Huabin Ge) Department of Mathematics, Beijing Jiaotong University, Beijing 100044, 
P.R. China 

E-mail: hbge@bjtu.edu.cn 


(Xu Xu) School of Mathematics and Statistics, Wuhan University, Wuhan 430072, PR 
China 

E-mail: xuxu2@whu.edu.cn 


12 



